What Is Claimed Is; 



1 LA method for using a computer system to solve a system of 

2 nonlinear equations specified by a vector function, f, wherein f(x) = 0 represents a 

3 set of nonlinear equations, / y (x) = 0,f 2 (x) = 0,f 3 (x) = 0, . . .,/„(x) = 0 9 wherein x 

4 is a vector (x I} x 2t x 3 , ... x„), the method comprising: 

5 receiving a representation of a subbox X = (Xj, X 2 , X„), wherein for 

6 each dimension, /, the representation of X l includes a first floating-point number, 

7 a h representing the left endpoint of X l , and a second floating-point number, b h 

8 representing the right endpoint of X x ; 

9 storing the representation in a computer memory; 

10 applying term consistency to the set of nonlinear equations, //(x) = 0, 

1 1 f 2 {\) = 0,/i(x) = 0, .,^(x) = 0, over X, and excluding portions of X that violate 

1 2 any of these nonlinear equations; 

13 applying box consistency to the set of nonlinear equations over X, and 

14 excluding portions of X that violate any of the nonlinear equations; and 

1 5 performing an interval Newton step on X to produce a resulting subbox Y, 

16 wherein the point of expansion of the interval Newton step is a point x within X, 

1 7 and wherein performing the interval Newton step involves evaluating f(x) using 

1 8 interval arithmetic to produce an interval result f ! (x). 

1 2. The method of claim 1 , wherein performing the interval Newton 

2 step involves: 

3 computing J(x,X), wherein J(x,X) is the Jacobian of the function f 

4 evaluated as a function of x over the subbox X; and 

5 determining if J(x,X) is regular as a byproduct of solving for the subbox Y 

6 that contains values of y that satisfy M(x,X)(y-x) = r(x), where 
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1 M(x,X) = BJ(x ? X) ? r(x) = -Bf(x), and B is an approximate inverse of the center of 

2 J(x,X). 

1 3. The method of claim 2, further comprising: 

2 applying term consistency to the preconditioned set of nonlinear equations 

3 Bf(x) = 0 over the subbox X; and 

4 excluding portions of X that violate the preconditioned set of nonlinear 

5 equations. 

1 4. The method of claim 2, further comprising: 

2 applying box consistency to the preconditioned set of nonlinear equations 

3 Bf(x) = 0 over the subbox X; and 

4 excluding portions of X that violate the preconditioned set of nonlinear 

5 equations. 

1 5 . The method of claim 1 , wherein applying term consistency to the 

2 set of nonlinear equations involves: 

3 for each nonlinear equation /(x) = 0 in the system of equations f(x) = 0, 

4 symbolically manipulating f£x)=0 to solve for an invertible term, g(x )\ thereby 

5 producing a modified equation g(x ))= h (x), wherein g(x )) can be analytically 

6 inverted to produce an inverse function g'\y); 

1 substituting the subbox X into the modified equation to produce the 

8 equation g{X)) = /z(X); 

9 solving for X) = g\h(X)); and 

1 0 intersecting X) with the vector element X } to produce a new subbox X + ; 
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1 1 wherein the new subbox X + contains all solutions of the system of 

12 equations f(x) = 0 within the subbox X, and wherein the width of the new subbox 

13 X + is less than or equal to the width of the subbox X. 

1 6. The method of claim 1 , further comprising: 

2 evaluating a first termination condition, wherein the first termination 

3 condition is TRUE if, 

4 zero is contained within f ! (x) ? 

5 J(x,X) is regular, wherein J(x,X) is the Jacobian of the 

6 function f evaluated as a function of x over the subbox X, and 

7 the solution Y of M(x,X) (y-x) = r contains X; and 

Q 8 if the first termination condition is TRUE, terminating and recording X as 

Q 

jg 9 a final bound. 

O 1 7. The method of claim 6, wherein the method further comprises: 

g ~ 2 evaluating a second termination condition; 

3 wherein the second termination condition is TRUE if a function of the 

2 4 width of the subbox X is less than a pre-specified value, s x , and the width of the 

O 5 function f over the subbox X is less than a pre-specified value, s F ; and 

6 if the second termination condition is TRUE, terminating and recording X 

7 as a final bound. 

1 8. A computer-readable storage medium storing instructions that 

2 when executed by a computer cause the computer to perform a method for using a 

3 computer system to solve a system of nonlinear equations specified by a vector 

4 function, t wherein f(x) = 0 represents a set of nonlinear equations, //(x) = 0, 
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5 f 2 (x) = 0,f 3 (x) = 0, . . . 5 /»(x) = 0, wherein x is a vector (x 1} x 2) x 3) ... x n ), the 

6 method comprising: 

7 receiving a representation of a subbox X = {X h X 2) wherein for 

8 each dimension, /, the representation of X l includes a first floating-point number, 

9 a l9 representing the left endpoint of X x , and a second floating-point number, b t9 

1 0 representing the right endpoint of X x ; 

1 1 storing the representation in a computer memory; 

12 applying term consistency to the set of nonlinear equations, fj(x) = 0, 

1 3 f 2 (x) = 0, f 3 (x) = 0, . . . , f n (x) = 0, over X, and excluding portions of X that violate 

1 4 any of these nonlinear equations; 

1 5 applying box consistency to the set of nonlinear equations over X, and 

1 6 excluding portions of X that violate any of the nonlinear equations; and 

1 7 performing an interval Newton step on X to produce a resulting subbox Y, 

1 8 wherein the point of expansion of the interval Newton step is a point x within X, 

19 and wherein performing the interval Newton step involves evaluating f(x) using 

20 interval arithmetic to produce an interval result f \x). 

1 9. The computer-readable storage medium of claim 8, wherein 

2 performing the interval Newton step involves: 

3 computing J(x,X), wherein J(x,X) is the Jacobian of the function f 

4 evaluated as a function of x over the subbox X; and 

5 determining if J(x,X) is regular as a byproduct of solving for the subbox Y 

6 that contains values of y that satisfy M(x,X)(y-x) = r(x), where 

7 M(x,X) = BJ(x,X), r(x) - -Bf(x), and B is an approximate inverse of the center of 

8 J(x,X). 
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1 10. The computer-readable storage medium of claim 9, wherein the 

2 method further comprises: 

3 applying term consistency to the preconditioned set of nonlinear equations 

4 Bf(x) = 0 over the subbox X; and 

5 excluding portions of X that violate the preconditioned set of nonlinear 

6 equations. 

1 11. The computer-readable storage medium of claim 9, wherein the 

2 method further comprises: 

3 applying box consistency to the preconditioned set of nonlinear equations 

4 Bf(x) = 0 over the subbox X; and 

5 excluding portions of X that violate the preconditioned set of nonlinear 

6 equations. 

1 12. The computer-readable storage medium of claim 8, wherein 

2 applying term consistency to the set of nonlinear equations involves: 

3 for each nonlinear equation /(x) = 0 in the system of equations f(x) = 0, 

4 symbolically manipulating /(x)=0 to solve for an invertible term, g(x thereby 

5 producing a modified equation g(x ))= h (x), wherein g(x can be analytically 

6 inverted to produce an inverse function g 7 (y); 

7 substituting the subbox X into the modified equation to produce the 

8 equation g(X)) =h(X); 

9 solving for X) = g \h(X)); and 

1 0 intersecting X) with the vector element X } to produce a new subbox X + ; 

1 1 wherein the new subbox X + contains all solutions of the system of 

1 2 equations f(x) = 0 within the subbox X, and wherein the width of the new subbox 

13 X + is less than or equal to the width of the subbox X. 

29 

Attorney Docket No. SUN-P7075-SPL Inventors: Walster et al. 

ARPH \SUN MICROSYSTEMS\SUN-P7075-SPL\SUK-P7075-SPL APPLICATION DOC 



1 13. The computer-readable storage medium of claim 8, wherein the 

2 method further comprises: 

3 evaluating a first termination condition, wherein the first termination 

4 condition is TRUE if, 

5 zero is contained within f ] (x), 

6 J(x 3 X) is regular, wherein J(x,X) is the Jacobian of the 

7 function f evaluated as a function of x over the subbox X, and 

8 the solution Y of M(x,X) (y-x) = r contains X; and 

9 if the first termination condition is TRUE, terminating and recording X as 
10 a final bound. 

>« 1 14. The computer-readable storage medium of claim 13, wherein the 

> I 2 method further comprises: 

O 3 evaluating a second termination condition; 

_T ? 4 wherein the second termination condition is TRUE if a function of the 

j*f 5 width of the subbox X is less than a pre-specified value, Sx, and the width of the 

y 6 function f over the subbox X is less than a pre-specified value, s F ; and 

Q 7 if the second termination condition is TRUE, terminating and recording X 

nj 

8 as a final bound. 

1 1 5. An apparatus that solves a system of nonlinear equations specified 

2 by a vector function, f, wherein f(x) = 0 represents a set of nonlinear equations, 

3 /y(x) = 0,/Xx) = 0J 3 {\) = 0, ...,f n (x) = 0, wherein x is a vector (x h x 2t x 3 , ... x„), 

4 the apparatus comprising: 

5 a receiving mechanism that is configured to receive a representation of a 

6 subbox X = (Xj 9 X% X n ), wherein for each dimension, i, the representation of 
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7 Xi includes a first floating-point number, a h representing the left endpoint of X x , 

8 and a second floating-point number, b l9 representing the right endpoint of X x ; 

9 a computer memory for storing the representation; 

1 0 a term consistency mechanism that is configured to apply term consistency 

1 1 to the set of nonlinear equations, fj(x) = 0, f 2 (x) = 0,f 3 (x) = 0, . ,.,f n (x) = 0, over 

12 X, and to exclude portions of X that violate any of these nonlinear equations; 

13 a box consistency mechanism that is configured to apply box consistency 

14 to the set of nonlinear equations over X, and to exclude portions of X that violate 

1 5 any of the nonlinear equations; and 

1 6 an interval Newton mechanism that is configured to perform an interval 

17 Newton step on X to produce a resulting subbox Y, wherein the point of 

18 expansion of the interval Newton step is a point x within X, and wherein 

19 performing the interval Newton step involves evaluating f(x) using interval 

20 arithmetic to produce an interval result f r (x). 

1 1 6, The apparatus of claim 1 5, wherein the interval Newton 

2 mechanism is configured to: 

3 compute J(x,X) ? wherein J(x,X) is the Jacobian of the function f evaluated 

4 as a function of x over the subbox X; and 

5 determine if J(x,X) is regular as a byproduct of solving for the subbox Y 

6 that contain the values of y that satisfy M(x,X)(y-x) = r(x), where 

7 M(x,X) = BJ(x,X), r(x) = -Bf(x), and B is an approximate inverse of the center of 

8 J(x,X). 

1 17. The apparatus of claim 1 6, wherein the term consistency 

2 mechanism is configured to: 
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1 apply term consistency to the preconditioned set of nonlinear equations 

2 Bf(x) = 0 over the subbox X; and to 

3 exclude portions of X that violate the preconditioned set of nonlinear 

4 equations. 

1 18. The apparatus of claim 1 6, wherein the box consistency 

2 mechanism is configured to: 

3 a Pply box consistency to the preconditioned set of nonlinear equations 

4 Bf(x) = 0 over the subbox X; and to 

5 exclude portions of X that violate the preconditioned set of nonlinear 

6 equations. 

1 19. The apparatus of claim 1 5, wherein for each nonlinear equation 

2 /(x) = 0 in the system of equations f(x) = 0, the term consistency mechanism is 

3 configured to: 

4 symbolically manipulate fi{\)=0 to solve for an invertible term, g(x )% 

5 thereby producing a modified equation g(x ))= h (x), wherein g(x )) can be 

6 analytically inverted to produce an inverse function g'\y); 

7 substitute the subbox X into the modified equation to produce the equation 

8 g{X)) = h(X); 

9 solve for X) = g \h(X)); and to 

1 0 intersect X) with the vector element Xj to produce a new subbox X + ; 

1 1 wherein the new subbox X + contains all solutions of the system of 

12 equations f(x) = 0 within the subbox X, and wherein the width of the new subbox 

13 X + is less than or equal to the width of the subbox X. 
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1 20. The apparatus of claim 1 5, further comprising a termination 

2 mechanism that is configured to: 

3 evaluate a first termination condition, wherein the first termination 

4 condition is TRUE if, 

5 zero is contained within f [ (x) 5 

6 J(x,X) is regular, wherein J(x,X) is the Jacobian of the 

7 function f evaluated as a function of x over the subbox X, and 

8 the solution Y of M(x,X) (y-x) = r contains X; and to 

9 terminate and record X as a final bound if the first termination condition is 
10 TRUE. 

i-JL 

0 1 21. The apparatus of claim 20, wherein the termination mechanism is 
'% 2 additionally configured to: 

1 jf 3 evaluate a second termination condition; 

Q 4 wherein the second termination condition is TRUE if a function of the 

g 5 width of the subbox X is less than a pre-specified value, ex, and the width of the 

fi 

JT 6 function f over the subbox X is less than a pre-specified value, e F ; and to 
7 terminate and record X as a final bound if the second termination 

O 8 condition is TRUE. 
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